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Cosmic Ray Proton Background Could Explain ATIC Electron Excess 
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We show that the excess in the Galactic electron flux recently published by Chang, et al. (Nature, 
20 Nov. 2008) can have a simple methodical origin due to a contribution from misidentifled proton 
induced electron-like events in the ATIC detector. A subtraction of the cosmic ray proton component 
from the published ATIC electron flux eliminates this excess in the range of 300 to 800 GeV. 



PACS numbers: 95.30.Cq, 14.80.Ly 

An excess in the Galactic electron flux in the range 
of 300 - 800 GeV recently published by J. Chang^ et al, 
[H has led to numerous speculations [l|, 0, 0, 0, IE B 
about the origin of the purported excess such as anni- 
hilation of dark matter (DM) P, |3k decaying DM [3,[3|, 
decay of lightest superparticle DM [j| , interaction of high 
energy cosmic rays with photon background near accel- 
erating sites 0] , a few nearby SNR Q , and distantyoung 
pulsars All theoretical predictions [3, S H, S 0, 01 
describe with more or less efficiency the published energy 
spectrum reported by Chang et al.,[l| in the energy re- 
gion of 300 — 800 GeV depending on a number of applied 
free model parameters. 

Chang, et al.,[lj use an excess of 70 events, spread 
over an energy range of 500 GeV, to report a ^ 6(7 sig- 
nal based on a 210-event sample, where a background of 
140 events was estimated with a GALPROP [8] calcula- 
tion, neglecting both systematic uncertainties associated 
with the GALPROP prediction jsj as well as those of 
the ATIC instrument. Not noted in reference [H, the 
expected fluctuation of the Galactic electron flux due to 
the stochastic nature of sources [3, [13] alone could intro- 
duce more than a 20% uncertainty for energies above 300 
GeV. A major source of background of the ATIC electron 
experimental data is due to the primary cosmic ray pro- 
ton component where the bulk of the proton energy is 
transferred to a single particle, for example, a 7r° in the 
carbon-scintillator target of the ATIC detector. 

In this letter, following the principle of Occam's Razor 
mentioned in the title of a recent paper by Stefano Pro- 
fumo , we investigated the methodical origin of the re- 
ported excess. We studied the contributions due to high 
energy primary protons and atmospheric and 7 to the 
reported Galactic electron flux in the ATIC data. The 
primary proton and atmospheric and 7-background 
components at the level of ATIC flights 5g/cm^) jlj 
and energies more than 50 GeV have an energy spectral 
power index 7^ — —2.7 and could change significantly 
the observed spectral index 7e ~ —3.2 of Galactic elec- 
trons, especially above the TeV energy region, where an 
exponential energy cut-off is expected [l3| ■ 

ATIC, the acronym for Advanced Thin Ionization 
Calorimeter was designed for multiple, long duration bal- 
loon flights, to measure nucleonic cosmic ray spectra 



from 10 GeV to near 100 TeV total energy, using a fully 
active Bismuth Germanate (BGO) calorimeter. It was 
equipped with the first large area mosaic of small fully 
depleted silicon detector wafers for charge identification 
of cosmic rays from H to Fe. As a particle tracking sys- 
tem, three projective layers of x-y scintillator hodoscopes 
were target [ll|. The experiment had no magnetic field 
and was designed for cosmic ray elemental composition 
measurements with limited particle ID (PID) capability 
especially for the types of measurements mentioned in 
reference [l[. 

In this letter, we used the published data of reference 
[H and applied their cut parameters, without delving 
into a detailed analysis of their validity, to a GEANT 
Monte Carlo (MC) [12^. The primary energy nuclei spec- 
tra were taken frorn power law approximations of balloon 
and satellite data llJI 
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where $a and 7a are the scale parameters and spectral 
indices for A = p, He, O-like, and i^e-like nuclear species 

The flux of the p, and 7-background components 
of cosmic rays at the ATIC level were computed using 
CORSIKA shower simulation code [15j in the frames 
of SIBYLL interaction model, taking into account 
the South Pole atmosphere and magnetic field. The es- 
timated contributions of the secondary atmospheric 7 
and e^-components to the Galactic electron flux varied 
slightly from the corresponding ATIC results of refer- 
ence PI, where the proton background was said to be 
included. This result has indicated that proton-induced 
c-likc events were substantially underestimated in refer- 
ence [l[. Because the contribution of and 7 atmo- 
spheric components are negligible to the overall data, we 
have omitted the details here and instead concentrate on 
the main background which comes from cosmic ray pro- 
tons. 

The contribution of the primary proton background to 
the Galactic electron flux was calculated using the proton 
spectra at the ATIC level in a GEANT-3.21 calculation 
with the GCALOR routines in the frame of the FLUKA 
interaction model for the ATIC detector [l3|. This cal- 
culation takes into account the ATIC geometry with the 
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BGO layer number o^+o^+fj+Fg (mm+mm ) 



FIG. 1: GEANT simulated electron (left panel) and proton 
(right panel) for 100 cascade events showing the longitudinal 
profiles in the ATIC BGO calorimeter at energies of 400 — 600 
GeV (upper lines) and at 50 — 100 GeV (lower lines). 



silicon matrix, scintillator hodoscopes, graphite targets 
and eight layers of x-y BGO crystals, containing appro- 
priate material compositions The total number of 
simulated proton and electron events were 100,000 and 
10,000, respectively. 

These GCALOR calculations for the ATIC detector 
were previously compared with a FLUKA-2006 calcula- 
tion 



17l | as well as with those of an ATIC test-beam 
run at CERN [i3|. Furthermore, we also ran FLUKA- 
2006 for the ATIC configuration and results were consis- 
tent with those of the GCALOR calculations. Figure [1] 
shows the electron (left panel) and proton (right panel) 
induced cascade profiles in the ATIC BGO calorime- 
ter after e-like cascade selection criteria: Wmax £ 5, 

Ed,ma^ > Ed,max+1 > Ed,maK+2 ■ ■ ■ Ecl,S: Ed,l > 2 GcV 

and Ed.s > 25 MeV, where Ed.n is the detected energy 
in the n-th BGO layers and -Ed, max — niax (i?rf_„) for 
n — 1, ... 8. The lower and upper lines in Figure [T] cor- 
respond to 50 - 100 GeV and 400 - 600 GeV detected 
energies respectively. The selection criteria rejects 74% 
of proton events with visible energies Ey > 50 GeV with- 
out significant change of the electron cascade detection 
efficiency. This result is in close agreement with the e- 
cascade selection efficiency (~ 90%) and proton rejection 
efficiency (~ 75%) of reference [ij. 

We simulated the lateral distribution of the cascade 
parameter F = ai + a2 + Ft + employed in ref- 



FIG. 2; Distribution of combined cascade parameter ai+a2 + 
F\ + F2 for detected energies Ebgo = Ea > 50 GeV (and 300 
GeV, inset histograms) . The symbols are the ATIC data from 
The large (line) and shaded histograms are the expected 
distributions for primary proton and electron passing through 
the BGO calorimeter. 



erence [ij for proton induced e-like cascades and pure 
e-cascades. Figure [2] shows the published data of ref- 
erence jl] and our calculations. Note the indices of a 
and F refer to the BGO layers. Our simulated data 
for Ey > 50 GeV was normalized to the correspond- 
ing statistics of the ATIC data fH taking into account 
energy-dependent ratio {d'^p/ dEp) / {d'^d dE^) of the ex- 
pected proton flux from expression (1) and the expected 
electron flux d^e/dE^ = aE~^ °'^ exp(— iJe/i?o) from the 
HESS experiment The results of the GEANT sim- 
ulations for visible energies Ey > 300 GeV are shown in 
the inset histogram, as well. 

It is observed that the contribution of primary pro- 
ton induced e-like cascades in the ATIC BGO calorime- 
ter to the flux of Galactic electrons is quite substantial 
in the region of the reported excess {Ed > 300 GeV). 

"Pure" e-cascade pattern recognition in the ATIC 
data is performed using diagonal energy dependent cuts 
[l[ in a 2-dimensional (2d) (cti + a2, Fi + F2) space. 
The corresponding GEANT simulated {{ai+a2), Ed) and 
{{Fj + Fs),Ed) scatter plots are presented in Figure [3l 
The number of simulated events was taken to be equal at 
5000 events for both electrons and protons with Ed > 50 
GeV. Note that the actual ratio of p/e in cosmic rays is 
approximately 300-500 for these energies. 
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FIG. 3: 2-dimensional ((cti + 0-2), Sbgo), (upper panel) and 
[[Fj + Fs),Ebgo), (lower panel) scatter plots for 5000 simu- 
lated proton and electron events in the ATIC BGO calorime- 
ter. The lines are the corresponding proton rejecting cuts. 



Because the 2d selection criteria (cuts) applied in the 
published data were unavailable from the literature, we 
derived these 2d cuts (lines in Figure[3]) from the GEANT 
simulated database using known 84% efficiency of e- 
cascade selection from [l[ . Lines shown in Figure [3] rep- 
resent cuts corresponding to the same afore-mentioned 
selection efficiency of e-cascade from GEANT simulated 
database. 

The energy spectrum of proton induced and subse- 
quent misidentified e-like events expected in the ATIC 
data were calculated according to the expression 



AiV 



Ntot AE-Se-^ 
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where Ee = Ed/s and e is the average fraction of re- 
leased energydue to e-cascade passing through the BGO 
calorimeter [l| , 3p is the integral energy spectrum of the 
primary protons from expression (1) for iSp^min — 100 
GeV, Ntot is the total number of simulated proton events, 
AiV and AE are the number of selected e-like proton in- 
duced events in the given AE energy bin, Se — 0.84 [l| 
is the efficiency of 2-dimensional e-cascade selection cuts, 
and ^(Ee) C = ^e,ATic /^e,Prim - 0.78 is the atmo- 
spheric reduction factor taken from reference [lO]. The 
corresponding results are presented in Figure H] (black 
circles) in comparison with published data of reference 
[H (black star symbols). The accuracy of the obtained 



energy spectrum was improved using the derivative of 
the corresponding integral spectra due to a priori known 
spectral index {-jp = —2.75) of proton induced e-like cas- 
cades. The corresponding differential energy spectrum, 
normalized for energy E^ = 100 GeV, is shown in Fig- 
ure m (black line with blue-shaded area for the expected 
errors), as well. Red star symbols are the residual en- 
ergy spectrum of Galactic electrons after subtraction of 
the expected proton induced e-like cascade from the re- 
ported data of reference [ij . 

Estimated contribution of the expected e-like proton 
background is turned out to be 96 ± 9.8 as compared 
to the reported 70-event excess of reference [H, thus re- 
jecting the conclusion of exotic electrons. 
In summary. Figure 2] clearly shows that the MC pro- 
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FIG. 4: Galactic electron and misidentified e-like energy spec- 
tra. The solid black stars (Fe) are the data of reference [l|], the 
red stars are the ATIC data after subtraction (Fe — fp-e) of 
e-like proton energy spectrum {fp-e, the solid line with blue 
shaded area of errors) obtained from the corresponding inte- 
gral flux, the solid black circles are the GEANT simulated 
proton induced e-like events energy spectrum. The dashed 
line is the power law energy spectrum with exponential cut- 
off from HESS [H. 

ton spectrum has a substantial contribution in the en- 
ergy region of more that 300 GeV where an excess of e- 
like events was reported [l[. Our analysis indicates that 
a calorimeter such as ATIC with its limited PID capa- 
bility could easily misidentify protons for electrons and 
thereby yield unwarranted physics conclusions. The esti- 
mate background of 96±9.8 events in the energy region of 
300 - 800 GeV completely explains the 70-event excess of 
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reference [l[ . We therefore conclude that the contribution 
from cosmic ray protons at energies above 300 GeV has 
not been properly included in reference [![. This analy- 
sis shows that the results of reference [H does not merit 
any exotic conclusion and that the electron spectrum is 
consistent with that reported by the AMS collaboration 
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